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We return to the question of how the choice of stabihzer generators affects the preservation 
of information on structures whose degenerate ground state encodes a classical redundancy code. 
ControUed-not gates are used to transform the stabilizer Hamiltonian into a Hamiltonian consisting 
of uncoupled single spins and/or pairs of spins. This transformation allows us to obtain an analytical 
partition function and derive closed form equations for the relative magnetization and susceptibility. 
These equations are in agreement with the numerical results presented in [Phys. Rev. A 80, 042313 
(2009)] for finite size systems. Analytical solutions show that there is no finite critical temperature, 
Tc— 0, for all of the memory structures in the thermodynamic limit. This is in contrast to the 
previously predicted finite critical temperatures based on extrapolation. The mismatch is a result 
of the infinite system being a poor approximation even for astronomically large finite size systems, 
where spontaneous magnetization still arises below an apparent finite critical temperature. We 
extend our analysis to the canonical stabilizer Hamiltonian. Interestingly, Hamiltonians with two- 
body interactions have a higher apparent critical temperature than the many-body Hamiltonian. 
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I. INTRODUCTION 

The equivalent of a magnetic memory for quantum 
information would consist of a macroscopic number of 
qubits with multi-qubit interactions that create a sin- 
gle stable qubit memory. The free energy of the system 
would depend upon an external control to spontaneously 
break global symmetry in the presence of environment- 
induced fluctuations. Kitaev's toric code in a four di- 
mensional lattice would achieve this task [l], 0] , but its 
implementation seems currently unlikely. Bravyi and 
Terhal have recently shown that a two-dimensional self- 
correcting quantum memory may not exist [Sj] . If dimen- 
sionality is an engineering limitation, the solution may 
be self-correcting memories of finite size based on con- 
catenated codes in which the number of qubits involved 
in each interaction grows with the lattice size [1]. The 
classical concatenated triple modular redundancy code in 
the formalism of quantum stabilizers using the standard 
choice of generators fulfills this prerequisite for classical 
memory. 

The stabilizer for a subspace is defined as the group 
of Pauli operators that act trivially on a code space and 
whose eigenvalues are -1-1. The code space is the de- 
generate ground state of a Hamiltonian built from the 
stabilizer elements with negative couplings. The triple- 
modular redundancy code is a textbook example for in- 
troducing the idea of stabilizer error correcting codes Q . 
Classical error correcting codes represent a subset of 
quantum error correcting codes that only protect against 
classical bit- flip errors but not phase errors @. At each 
level of concatenation k, the logical bit consists of three 
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bits of level fc — 1, and correction works by majority vote 
at the lowest level first and then working up. The fc-th 
level of concatenated code contains 3*^ bits or classical 
spins, and it can always correct a maximum of 2*^ — 1 
errors on the physical bits. The increase of k leads to 
many-body operators that test the parity of | x 3'^ bits 
at once. This exponential increase in the many-body na- 
ture of the Hamiltonian makes the physical construction 
of such a system unrealistic. 

An alternative choice uses only elements that test a 
pair wise agreement. This set of Pauli operators gener- 
ates the same stabilizer group and represent an Ising 
Hamiltonian with characteristic thermodynamic and ki- 
netic properties. Using Monte Carlo simulations, we ex- 
amined the thermal magnetization of this pairwise choice 
of stabilizers (Structure 1 in Fig. [T]) and the effect of 
adding non-independent stabilizers to the Hamiltonian 
(Structures 2 and 3) [6j. For Structure 1, 3*^ — 1 inde- 
pendent stabilizer elements form a tree. Structures 2 and 
3 are modifications that include cycles in the structure. 
The cycles are equivalent to choosing an overcomplete set 
of stabilizer generators. 

In this paper, we analytically evaluate the choice of 
stabilizer generators on the preservation of information 
. Specifically, a unitary operator is constructed from 
controlled-not gates that converts a Hamiltonian repre- 
senting an Ising tree into a Hamiltonian of uncoupled 
spins in a magnetic field. Applying the same unitary 
operator to the tree-like graphs of Structures 2, 3, and 4 
yields partition functions corresponding to a collection of 
independent single spins and independent pairs of spins. 
A slight modification of the sequence allows us to calcu- 
late the analytical magnetization of the canonical stabi- 
lizer Hamiltonian. The results presented here agree with 
our previous numerical work for relatively small, finite- 
size systems. Closed form partition functions for each 
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size scaling fails to describe magnetic susceptibility peaks 
shifted away from Tc—0. For Ising trees and for Sierpin- 
ski gaskets, the relative magnetization approaches zero 
in the thermodynamic limit 0, [l^ - [l6j , but it persists for 
very large systems (comparable to the number of hadrons 
in the universe) [1, [TtI . The nature of the magnetic phase 
transition for an infinite system is not applicable to sys- 
tems of laboratory dimensions. We find similar behav- 
ior in Structures 2, 3, and 4, but with higher apparent 
critical temperatures (defined as the temperature where 
magnetic susceptibility reaches its maximum). Surpris- 
ingly, the canonical choice of elements to generate the 
concatenated three-bit error-correction code exhibits the 
lowest of the finite-size apparent critical temperatures. 
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FIG. 1. Transformation of the memory stabilizer structures 
generated by two body interactions from the Ising basis to 
the free-spin basis. Black dots and open circles are spin sites 
(qubits), and the lines show pairs of interacting spins (gener- 
ators). In the free-spin basis, the black dot without interac- 
tions is a single free spin, open circles are independent spins 
in a magnetic field, and the connected circles are independent 
pairs of interacting spins in a magnetic field. The interaction 
strength J is constant (see Eq. [T|. The total number of bits 
increases with concatenation level, k, as 3''. Only fc = 3 level 
structures are shown. 



of the four self-correcting memory structures allow us to 
examine the problem at much larger k. 

A direct measurement of the degree of preservation of 
the information can be read from the spontaneous mag- 
netization at zero magnetic field. Below a certain temper- 
ature, a single spin, sq, is sufficient to bias the system into 
one of the two states of broken symmetry. The finite-size 
system develops spontaneous magnetization and the sin- 



gle order parameter mp 



(Sjlo^ {soSj))/N approaches 



the value of 1 ■ The stability of the structure, as mea- 
sured by the temperature range in which mo is preserved, 
depends on the energy barrier that separates the two 
ground states and the number of pathways that traverse 
the barrier. 

Structure 1 is an example of an Ising tree with free 
boundaries. The Ising model on Cayley trees results in 
partition functions that are equivalent to free spins Q . 
Our previous analysis, based on TV — 81,243,729, and 
2187 bits (fc = 4 — 7) and under the assumption that 
Fisher's finite-size scaling method [l3| applies to these 
type of Ising graphs, yielded a non-zero Tj,. But contrary 
to Sierpinski fractals [ll| - [l3| . where a few data points 
seem to be enough to forecast Tc correctly, the finite- 



II. CNOT TRANSFORMATIONS AND ISING 
SYSTEMS 



The algebra of controlled-nots (CNOTs) and Pauli Z 
operators from quantum computation is used to find 
analytical solutions for the internal energy and mag- 
netization of the Ising structures in Fig. [1] Follow- 
ing standard notation, the spin or qubit basis is labeled 
|0) and |1) with the Pauli Z operator in the computa- 
tional basis acting as Z \x) — (—1)^ \x), where x equals 
or 1. The controlled-not operation on two qubits 
can be written compactly in the computational basis as 
CA^Or(l,2) |a;2) = \xi) \x2 ® xi) where qubit 1 is 
the control qubit and © represents addition modulo 2. 
Through out this manuscript, we take advantage of the 
following relations: 



CNOT{j, k)CNOT{i, k)=I 
Z,Z, = I 



CNOT{j, k)ZjCNOT{j, k) 



CNOTij,k)ZkCNOT{j,k) 
CNOT{j,k)ZjZkCNOT{j,k) = Zu 



ZjZk 



The last two relationships convert between Ising cou- 
plings, ZjZ}^^ and local magnetic fields, Z^. The re- 
peated application of CNOT transformations is an ex- 
plicit method to obtain the zero-field partition function 
for any Ising tree Hamiltonian of N spins, which is always 
equivalent to a single free spin and — 1 independent 
spins in a magnetic field [9]. The same transformation 
applied to trees that are graphs with a few cycles re- 
sults in partition functions of clusters of spins. All of 
these partition functions are products of partition func- 
tions of few spins and do not lead to any singularities of 
the zero-field thermodynamic response functions. They 
do, however, lead to differences in magnetic behavior. 
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III. ANALYTICAL SOLUTION FOR THE 
PARTITION FUNCTION 

The stabilizer is defined as all the products of Pauli 
operators that act trivially on the code space. For the 
bit-flip code, we can choose any set of pairwise Ising in- 
teractions that generates the stabilizer operators. These 
generators form a Hamiltonian that is similar to the fer- 
romagnetic Ising model, 

H=-jY^Za,, (1) 

(id) 

where indicates a sum over nearest neighbors, and 
J sets the energy scale of the problem with temperature 
measured in units of J/Zcs. The choice of generators 
determines the structure and properties of the system 

a- 

1. Ising trees 

A tree is a connected graph without cycles or loops. 
As a result, there is one and only one path between any 
two nodes. In an Ising tree, the nodes represent bits and 
the edges represent the Ising interaction. Each node, 
n, is connected to a single parent, Up, and one or more 
children, n,.. If the node n is at a distance d from the 
root, the parent is at a distance d—1, and the children are 
at a distance c? + 1. For convenience, we define a function 
D that converts labels to the minimum distances from 
the root, e.g., if D{ri) = d then D{nc) = d-\-l. We label 
each node by its number and its parent's number to make 
explicit the tree nature of the graph. The Hamiltonian 
for TV spins is then written as 

N-l 
n— lie 

and the Z operator on the root is labeled ^[o,o] although 
the root has no parent. 

The CNOT{[np,n],[n,nc\) operator transforms 
Z[n^,n\Z[n,ne] iuto Z^^^ri,], but also transforms 
^[n,n,]^[n,,ng,] into Z[„p _„] Z[„_„^] Z[„^_„g^] , where gc 
labels the children of the children. By applying CNOTS 
first at the outermost connections (leaves) and then 
moving inwards, we can effectively transform all of the 
Ising terms into single spin terms. 

We define 

^= n n CNOT{[np,n],[n,n,]) (3) 

d=0 D{n)=d 

and the product implies right multiplication. Applying 
this unitary to the Hamiltonian of Eq. [2] yields 

H' = UHU'' 

N 
n=l 



which represents N — 1 spins in a magnetic field and one 
free spin. We will refer to this basis as the free-spin basis 
and the original computational basis as the Ising basis. 
U is the transformation matrix between the two bases 
(see Fig. [J). 

In Ising trees, every qubit, except the root, has 
the Hamiltonian Hi — —JZ in the free-spin basis. 
The partition function is then simply the product of 
the partition function of a single spin in a magnetic 
field: Qi = exp ( J/ftsT) exp (- J/fcsT). The ther- 
modynamic density matrix for a single spin is p = 
l/2[/2x2 + tanh( J//cbT)Z]. The density matrix is used 
to calculate the polarization in the free-spin basis, e = 
Tr [Zp] = tanh( J/fcs T"), and the internal energy, < 
Ex -JTrfZp] = -Je. The total thermal density 
matrix for all of the N spins is a tensor product over 
independent spin density matrices, 

Ptotal = ®n=0Pln^,n] = ^2x2 <^n=l (5) 

and the total internal energy is then < Etotai >= 
Tr[H' Ptotal] = -J{N - l)e. 



2. Ising trees with cycles 

The CNOT Ising tree transformation can also be ap- 
plied to graphs that can be decomposed into a spanning 
tree and Ising couplings between spins with the same 
parent (siblings). The Hamiltonian for N spins is now 

N-l 

^ = ^"^ X/ X/ ^["p.n]-^[n,"c] ^ J X/ ^[np,n]-^K,m] 
n=0 Me (n.m) 

(6) 

and the same unitary of Eq. [3] transforms it to the free- 
spin basis, thus 

N 

H' ^ Z[np,n] - J ^ ^K,n]2'[„p,„i]. (7) 

n=l (n,m) 

This is the Hamiltonian of one free spin and finite Ising 
graphs of sibling spins in non-zero magnetic field. 

Here we examine connections only between sibling 
pairs, that is the triangular cycles in Structures 2, 3, 
and 4 (Fig. [1]). In this case, there are three types of 
spins: i) the root which is depolarized in the free-spin 
basis and has {Eq) = 0, ii) the spin that is not connected 
to a sibling and is described by Hi ~ —JZ, which is 
equivalent to a spin in a magnetic field, and iii) spins 
that are connected to a sibling that have the two-spin 
Hamiltonian H2 ~ —J{Zi+Zj + ZiZj). The expected 
energy of the spins with the magnetic field Hamiltonian is 
(El) = —Je with magnetization e — tanh{J/kBT). The 
partition function of the siblings is Q2 = exp {3J/kBT) + 
3exp (—J/kBT), and the two-spin density matrix is then 
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l/4(/42;4 + oiZi + aZj + aZiZj), where 

_ exp {■jJ/kBT) - exp {-J/kBT) 
" ~ exp (SJ/ksT) + 3 exp {-J/ksT) ' 



(8) 



The energy is < E2 >= — 3 Ja, and the magnetization of 
a single spin is a. 

The total internal energy is the sum of energies for 
the three types of spin, < Etotai >=< Ei > Ni+ < 
E2 > N2/2 = -J{eNi + faiVs). The internal energies 
for Structure 1, Structure 2, Structure 3, and Structure 
4 are then -Je{3'' - 1), -J[e{3''-^ - 1) + |a(2 • 3'=-^)], 



-J|a{3'= - 1), and -J[e(2 • 3''-^) + |a(3'=-i - 1)], re- 
spectively. 

The expectation value of the operators constructed 
from products of Z's can be calculated quickly from the 
density matrices for the three spin types. The root is 
unpolarized, po — \/2l2x2, and as a consequence any op- 
erator that contains 2^[o,o] '^iH be zero. The single spins 
will contribute e per Z . The paired spins are correlated 
and will contribute a for individual Z's {Zi, Zj) or the 
product (ZiZj). These rules are sufficient to calculate 
the magnetic properties of the system and have a suc- 
cinct description in terms of the geometry. 



A. Calculation of the magnetization and the 
magnetic susceptibility 

The magnetization operator in the computational or 
Ising basis is M = '^^Zq Z^^p^n] and its expectation 
value is zero by symmetry. The product of the mag- 
netization of each spin and the magnetization of the root 
define the relative magnetization operator 



N-l 



(9) 



The root spin is sufficient to bias the system into one 
of the two states that break the symmetry Q, and the 



Mj /N is non-zero in the thermodynamic limit when 

the system is in a ferromagnetic phase. The square of 
the magnetization relates to the magnetic susceptibility 
per spin as follows: 



M 



My 



NkgT 



NkeT 



(10) 



Each of the Zn^^n operators must be transformed into 
the free-spin basis in order to calculate the magnetic 
properties. The basis transformation of Eq. [3]maps each 
Znp,n operator onto a product of Z's. When n is at a 
distance d from the root, the transformation yields 

= ^[nrf_i.n]^[«rf_2,nd-ll---^[0.ni]^[0,0], (H) 

with each parent labeled as n^^i. The operator Zj^^ „] 
becomes a product of Z's on every node on the path from 
the root to the spin n. 



(a) ® (b) ® 



a a 
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FIG. 2. Example of labeling edges in (a) trees and (b) tree-like 
graphs to calculate the magnetic properties based on paths 
between nodes (see text). 



The key observation is that the local magnetization 
operators in the Ising basis are transformed into paths in 
the free-spin basis. Calculations can then be performed 
using the paths as follows: 

• Label the edges of the tree-like graph with paired 
siblings by a, if the edge is part of a triangle, or by 
e, otherwise. 

• Define Path(n, I) as the product of the edge labels 
between nodes n and I along the shortest path. 

Fig. [5] shows a tree and related tree-like graph 
with the edges labeled. As an example we calcu- 
late Path(4,9). For the tree (Fig. EK), Path(4,9) = 
Path(4, l)Path(l,5)Path(5,9) = e^. In the tree- 
like graph (Fig. [5]d), there is a shortcut between 
the paired sibling nodes 4 and 5 and Path(4, 9) = 
Path(4, 5)Path(5, 9) = ea. 

As shown in Appendix |Al the magnetic thermody- 
namic averages can be related to the paths as 



N-l 



(M) = Path(0, 



and 



N~l N-l 



n=0 1=0 

For trees, these expressions simplify to 



(M) = J2 fidy 



and 



d=0 



2d„ 



(12) 



(13) 



(14) 



(15) 



where f{d) is the number of nodes a distance d from the 
root and 4>{d) is the number of unidirectional paths of 
length d. 

In summary, notice that in the Ising basis, H encodes 
the geometry by selecting which spins are paired (Eq. [TJ , 
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and that the magnetization operator is independent of 
the connectivity of the N spins. In the free-spin basis, 
H' is independent of the graph for trees with N nodes 
(Eq. 2]), and the geometry is now encoded in the magne- 
tization operator (Eq. [12]). 

This is weU illustrated by calculating the magnetiza- 
tion for two simple examples: a line of iV-spins and iV — 1 
spins connected to a central spin. In both cases the trans- 
formation to the free-spin basis results in a Hamiltonian 
of TV — 1 spins in a magnetic field and a single free spin. 
As a result, the partition function and density matrix in 
the free-spin basis are equivalent; however, the magneti- 
zations are quite different. For a line, /(d) = 1, therefore 

^M^ — J2d=o ^'^ converges to 1/(1 — e) in the limit of 
large N . This yields the familiar result that the magne- 
tization per spin is vanishingly small for T > 0. For the 
central spin case, c? = 1 or and /(I) = iV — 1, with 

the resulting magnetization being {^M^ = [N — l)e + 1. 

The system has non-zero magnetization per spin for all 
T < oo. 

Below we use the equations derived in this section to 
find analytical expressions for the magnetization and the 
susceptibility of the stabilizer structures of Fig. [T] All of 
these systems grow in size as iV = 3*^ as they are based 
on the concatenation of three units of 3*^""'^ spins at each 
level k. The path from the root to the furthermost spin 
is of length dmax = k. 



One can understand the result by imagining building up 
the tree level-by-level. The level k adds 2 nodes to every 
node in a level k — 1 tree. The paths between nodes 
and the root in the inner k — 1 tree are the same, and 
the leaves add two paths that are one edge longer. This 
results in the recursion formula: 



M(k)) 

/SI 



(l + 2e)(M(A:-l) 



SI 



(18) 



Notice that this last equation also generates Eq. [TTl thus 
the relative magnetization per spin at zero magnetic field 
is 



nT-osi(fc) 



Msi{k) 



l + 2e 



(19) 



N V 3 

which vanishes in the limit of large k for all e < 1 and 
T > 0. In order to calculate the magnetic susceptibility 
using Eq. I10[ we need to first evaluate the magnetization 
squared operator. Starting from k — \, two leaves are 
added to every node. Each path of length d > on the 
k — \ tree now has two extra leaves on each end. This 
results in one path of length d, four paths of length d -f 1 , 
and four paths of length d + 2. For the 3*^"^ paths of 
d = 0, there are now two paths of length one, one path 
of length two, and two new paths of zero length. Using 
these observations and defining (d, k) as the number 
of paths of distance d between two spins, Eq. [T2]can be 
written recursively: 



1. Structure 1 

In the Ising tree labeled Structure 1 of size iV = 3'^, 
the number of nodes at distance d from the root is 



(16) 



and according to Eq. [TH the expected relative magneti- 
zation of Structure 1 at level k is then 



2k 



(M2(fc))^^ =3'^ + ^20i(d,A;) 



d=l 

i''-'^ (l + 2 + 4e + 2e2) 

2(k-l) 

-h(l + 4e + 4e2) ^ 2(/)i(d,fc-l) 

d=l 

{l + 2ef {M\k^ I)) 

+2(1 - e2)3fe-i. (20) 



M{k)j^ " X! ( ^ ) " + (-^'^^ ^^'^ solution to the recursion formula is 



{M^{k))g^ = {l + 2e)2'= + 2(1 - e^){l + 2e)2('=-i) 



l-[i/{l + 2eff 
l-3/(l + 2e)2 



(21) 



and the magnetic susceptibility per spin is then 



2. Structure 2 



Xsi{k) 



2(1 - e2)(i + 2e)2('=-i){l - [3/(1 + 2e)2]'=} 



7VfcBT[l-3/(l + 2e)2] 



(22) 



Structure 2 is similar to Structure 1 but the leaves 
are connected forming triangular cycles. The number of 
spins at the minimum distance d from the root is the same 
as in Structure 1 but now there are single spins and spin 
pairs in the free-spin basis. For paths that include leaf 
spins from Structure 1, the magnetization needs to in- 
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elude the polarization of a spin pair, a. Structure 2 with 
3*^ spins is equivalent to Structure 1 with 3*^^^ spins with 
a sibling pair connected to each spin. The magnetization 
is then 



M(fc))^_^ = {l + 2a){Msi{k-l) 



(23) 



The thermodynamic average of (Af^) for a Structure 
2 of S'^ nodes can be built from a Structure 1 with 3*^"^ 
nodes by examining the extra shortest paths due to the 
attached outer cycles. The main difference is that the 
two new nodes connected to the k — 1 structure are a 
distance 1 apart instead of a distance 2. The result is 
that 



{M^k))^^^{l + 2af{M\k-l))^^ 
+2-3'=^i (l + a-2a2) 



(24) 



and 



XS2{k) 



(1 + 2a)^xsi{k - 1) + 2 (1 + g - 2a^) /keT 
3 

(25) 



3. Structure 3 

In this structure, each spin is part of a triangular cy- 
cle, and all spins but the root are paired spins. The 
thermodynamic average of the magnetization is identical 
to Structure 1 except the polarization is now a instead 
of e. 



M(k)) ^(l + 2a)* 

S3 



(26) 



The magnetization squared depends on the number of 
shortest paths between all spins, which is quite different 
from Structure 1 due to shortcuts made by triangular 
cycles. Applying the same building method of adding 
nodes to the core yields the following recursion relation 

(M2(fc))^3 = (l + 2a)2(M2(fc-l))^3 

+2(l + a-2a^)3''~\ (27) 

whose solution is: 



(A/2(fc))^3 = (1 + 2a)2'= + 2(1 + a - 2a''){l + 2a) 



2(k-i) f l-[3/(l + 2a)2]fc 
l-3/(l + 2a)2 



(28) 



The magnetic susceptibility is then 

2(l + a-2a2)(l + 2a)2(fc-i) 



Xsaik) 



NkgT 

{l-[il{l + 2afn 
[l-3/(l + 2a)2] ■ 



(29) 



4- Structure 4 

In Structure 4, each of the spins form part of triangular 
cycles except the outer nodes. The relationship between 
Structure 4 and Structure 3 is similar to the relationship 
between Structure 2 and Structure 1, and the magneti- 
zation properties are calculated to be 



M{k) 



Si 



M(fc-l)) (14- 2e), (30) 

S3 



{M\k))^^^{l + 2ef{M\k-l))^^ 

+2 ■ i^'^l - , (31) 



and 



. (l + 2e)2xS3(fe-l) + 2(l-6^) IksT 
Xsi{k) = ^ '- .(32) 



B. Extension to the Canonical Stabilizers 

The stabilizer formalism of quantum computing defines 
a subspace of n-qubits by a set of commuting observ- 
ables that are products of Pauli matrices on the n-qubits 
and have the value of 1 on the subspace. The stabilizer 
generators are independent operators, trace orthogonal, 
and commute with one another. As a result, there is 
always a unitary transformation which maps the stabi- 
lizer elements to Z operators on independent spins. Fur- 
thermore, this unitary can be constructed from CNOTs, 
Hadamards, and Pauli matrices Q. 

Structure 1 is derived from the three-qubit classical 
stabilizer code. The choice of generators is chosen to 
form an Ising tree, and this is not the standard choice. 
The standard choice is to use logical Ising interactions at 
every level of encoding. This choice results in generators 
that are multi-qubit interactions which grow exponen- 
tially with the level of encoding. 

The transformation that takes the stabilizer elements 
to independent Z's is closely related to the transforma- 
tion used for Structures 1, 2, 3 and 4. Instead of simply 
applying the CNOTs with the control on the inner node 
and then progressing inward, the control is alternated 
from inner to outer. A comparison of the two transfor- 
mations is shown for 9 qubits in Fig. |21 For Structures 
1, 2, 3, and 4 only A and B are applied. For the full 
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FIG. 3. A description of U as a quantum computing circuit 
for 9 qubits. For the tree and tree-like structures examined, 
the CNOTs are applied with the control towards the root 
starting from the leaves and then moving down layers until 
the root (A,B). For the full-stabilizer, the direction of control 
is alternated before applying the CNOTs at the next layer (A, 
A', B, B'). 



stabilizer, A, A', B, and B' are all applied. The detailed 
description of the transformation can be found in the Ap- 
pendix [B] The expected relative magnetization and the 
magnetization squared are respectively: 



M{k)) ^l + 2e 



3]k 



l_[(2 + .)e3] 
l-(2 



and 



where C = (2 + e)e. 



1 - (CV3)^ 
1 - (CV3) 



(33) 



(34) 



IV. RESULTS AND DISCUSSION 

A. Apparent Critical Temperature for Finite Size 
Systems 

The analytical results obtained from Eqs. [221 US] and 
1291 matchperfectly with our previous Monte Carlo sim- 
ulations f6l|. As an example. Fig. |4] compares the closed 
form equation of the magnetic susceptibility for Struc- 
ture 3 with numerical simulations of systems of various 
sizes. The susceptibilities are calculated using the ther- 
modynamic statistics of 5 x 10"* independent spin con- 
figurations generated with Wolff cluster simulations at 
different temperatures. 




FIG. 4. Magnetic susceptibilities per spin as a function of 
temperature T in units of J/fcs for different concatenation 
levels of Structure 3. The solid lines are calculated from Eg. 1291 
and the symbols are the result of Monte Carlo simulations [y] 
using 5 X 10* independent spin configurations. 



For finite-size self-correcting memories, the susceptibil- 
ity as a function of temperature shows a maximum which 
occurs at an apparent critical temperature T^g^^{N). It 
is clear from Fig. 2] that this temperature decreases with 
system size as expected. Finite size effects replace the 
divergences at the thermodynamic critical point by fi- 
nite peaks shifted away from [l^, [13 • Previously @ , 
we used a first order approximation to estimate these 
shifts for the case of susceptibility. A fit of T^^^ against 
the system size N gave us an estimate for Tc, xo and 
v' . With only four data points, we forecasted a finite 
Tc- Analytical solutions for the magnetic susceptibility 
permit the study of bigger systems and present a more 
complete picture of the finite-size effects on the magnetic 
properties. We obtain T^^,^{N) numerically and Fig. [S] 
compares them for the four Ising stabilizer structures, 
the canonical stabilizer, and the ID Ising model as a 
function of total number of spins in a double log scale, 
logjQ(log3(N)). The numerical calculation of T,^^^ for the 
ID Ising model of systems bigger than 3^^ spins results 
in a numeric underflow. The dotted line in the figure is 
an extrapolation using a two parameters fit of the solid 
line to the equation T^^^= ak~^. 

In the limit of systems of infinite size, Eqs. [T71 [231 [HI 
I30[ and [33] reveal that the only temperature at which mg 
takes the exact value of one is Tc— 0. It is seen in Fig. [S] 
that T^^^ converges very slowly to zero with system size. 
Based on the closed form equations this prolonged decay 
cannot be captured by a simple first order equation in 
N~^^^ nor by any finite power expansion in N with- 
out including an offset. We cannot calculate numerically 
the size of memory stabilizers with a T^^^ of practically 
zero before we run into numerical overflow. As shown in 
the figure, memory stabilizers utilizing all the observable 
matter in the universe will still behave as a finite-size sys- 
tem with almost all of their spins correlated at a finite 
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FIG. 5. Temperature of maximum relative magnetic suscepti- 
bility T^ax for different memory stabilizers of sizes that span 
from tens to 10^^ spins. Note that 1.70 and 2.23 in the ab- 
scissa correspond respectively to the Avogadro's number and 
to the predicted number of hadrons in the observable uni- 



apparent critical temperature on the order of J/ks- This 
is in contrast to the linear spin case where T^^j, rapidly 
approaches zero with increasing system size. For 3 spins, 
Structure 1 and the line are equivalent with T^g^^= 1.07. 
We can then ask how many spins are required to reach a 
certain T^^^. As an example, a maximum susceptibility 
of T^^^— 0.29 is achieved using 3^ spins in a line, but it 
would require a Structure 1 of 3'^^'^ spins. This presents 
an interesting challenge as these networked spin systems 
stand in contrast to our standard notion of what size the 
thermodynamic limit is appropriate. 



The thermal stability of the information encoded 
into finite systems for the four structures and the full- 
stabilizer can be related to T^^^ (Fig. [5]). The choice of 
stabilizers leads to a significant change in this apparent 
critical temperature. Two thirds of the spins in Structure 
2 form closed cycles, and the other third of spins form a 
core that is the same as a fc — 1 Structure 1. As expected 
for small systems. Structure 2 remains magnetized for a 
broader range of temperatures than Structure 1. Simi- 
larly, Structure 4, with 2/3 of its spins as free leaves, is 
less stable to temperature driven fiuctuations than Struc- 
ture 3. The cores of Structures 1 and 2, and Structures 3 
and 4, account for 1/3 of the spins, and they show similar 
magnetic behavior in the astronomical limit of large k. 
The canonical choice of stabilizer elements remains mag- 
netized below a broad range of temperatures, and shows 
the same long ranged order properties, but under this 
thermodynamic criteria, is a less efficient memory stabi- 
lizer than the simpler pairwise interaction geometries. 



B. Power-law Correlations and Finite Size Effects 

Assuming that there is a relation between a 1-D path 
of correlated spins and the total size of the system, 
L = iV^/'*, we can define a correlation length exponent 
scaled to the system size v' — v ■ d. The dimension, d, of 
each of the structures of Fig. [1] is unknown. According to 
the standard scaling hypothesis, and provided that the 
system size is large enough, the following scaling prop- 
erties are expected at the critical point: c{N) oc iV~, 

m{N) oc N-^, and x{N) oc [S]- Closed form 

equations for the magnetization and the magnetic sus- 
ceptibility can be used to calculate critical exponents. 

An analytical formula for the relative magnetization 
critical exponent can be obtained by equating the rela- 
tive magnetization per spin to the N~'^ power law. Af- 
ter simplification of the exponent fc, the P /v' critical ex- 
ponent can be written in terms of the logarithm of the 
relative magnetization as follows: 



P/v' = 1 - 



ln(3) ' 



(35) 



where ln(?/') = ln((M))/fc. For Structures 1 and 3, ■0 is 
independent of k and we find ipi = l-|-2e and V'3 = 1 + 2q!, 
respectively. For Structures 2 and 4, -0 depends weakly 

on fc as 0-2 = (If)'^' and 04 = 03 The 
power law for the relative magnetization per spin can 
thus be written as a function of any arbitrary tempera- 
ture: 



mo(T) = iV- 



MjKr» 1 

ln(3) J 



(36) 



This last equation reduces the magnetic susceptibility of 
Eqs. [ll]and[ig]to the form: 



x(fc,T) = A(fc,T) • (3'^)-2/3/^'+i^ 
where for Structures 1 and 3 A is 

(l-e2)[(i + 2e)2fe-3'=] 



A{k,T)si 



2k 



and 



A{k,T)s3 



fcBr(2e2-f 2e-l)(l + 2e) 



(l + a~2a^)[(l + 2a)^fc-3''] 
kBT{2a^ + 2a - 1)(1 + 2a)2'= 



(37) 



(38) 



(39) 



respectively. As T approaches the thermodynamic crit- 
ical temperature of Tc = 0, the A function of Eqs. [38] 
and [35] is independent of the system size and converges 
respectively to 



AiT) 



SI 



(1-a 



/fcBT(2e2 + 2e- 1) 



and 



AiT) 



fl 



S3 



2a2) 



kBT{2a^ + 2a- 1) 



(40) 



(41) 
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FIG. 6. Power law proportionality function A{k, T) for (a) 
Structure 1 and (b) Structure 3. Solid lines are calculated 
using Eqs. [38] and [39] while the circles come from the approxi- 
mations of Eqs. l40l and Bn In a broad region of temperatures 
near the thermodynamic critical point, these simple equations 
are enough to define a magnetization power law on for any 
system size. 



In this limit, the magnetic susceptibility can be writ- 
ten as the well known formula x{^) — ^{T) ■ N~ . 
Examining Eq. [HTl we find that j/i^' -I- 2/3/i/' = 1, 
which matches the Rushbrooke and Josephson scaling law 
d = j/v + 2/3/1' a written as a function of the correlation 
size exponent z^'. Fig. [Sj tests the temperature region 
in which this approximation holds for all system sizes. 
There is a broad temperature region where the magnetic 
properties are well described by N and temperature de- 
pendent critical exponents. 

The set of critical exponents obtained numerically for 
particular ill-predicted critical temperatures (see Tables 
IV and V of Ref. [6] ) match very well with those calcu- 
lated using Eq. [3S] and the hyperscaling relation. These 
apparent critical temperatures fall in the temperature re- 
gion in which Eqs. HHl and UT] hold. 

We find that there is a broad temperature region above 
Tc where the /3/i>' exponent is almost zero (it is strictly 
zero only at T = 0). Complementary, and for the same 
broad temperature region, the j/v' critical exponent 
reaches almost the value of one. The interpretation is 
simple: the number of correlated spins grows almost at 
the same rate as N. Spins that present long ranged cor- 
relations develop a net macroscopic alignment when an 
infinitesimal magnetic field is applied [7|. 



V. CONCLUSION 

We use CNOT gates to transform the Hamiltonian of 
stabilizer structures into a Hamiltonian consisting of un- 
coupled single and pairs of spins. In the original basis, 
the Hamiltonian encodes the geometry by selecting which 
bits interact. The magnetization operator is independent 
of the graph for N spins. In the free spin basis, the Hamil- 
tonian is independent of the graph for spins and the 
geometry is now encoded in the magnetization operator. 
This transformation allows us to obtain an analytical par- 
tition function and closed form equations for the effective 
magnetization and susceptibility with respect to a central 
spin. 

The analytical solutions match very well with the nu- 
merical results presented previously 16| for finite size sys- 
tems of TV < 3® spins. With a slight modification of 
the transformation sequence we calculate the analytical 
magnetization and susceptibility of the canonical stabi- 
lizer Hamiltonian. 

In our previous calculation based on four values of fc, 
we forecast a finite critical temperature for systems of 
infinite size. Our analytical solution shows that this pre- 
diction is incorrect. After applying a sequence of CNOT 
operations on the four stabilizer Hamiltonians studied in 
this work, the partition function results in a collection of 
free elements. The interactions represented in the mag- 
netization operator yield to graphs that have a transition 
from magnetic to random at Tc=0 in the thermodynamic 
limit. However, they possess unusual long-range-order 
properties as previously observed in hierarchical systems 
which also have Tc=0 (e.g.: Sierpinski gaskets |T7| and 
Cay ley trees Q). The memory stabilizer structures de- 
velop spontaneous magnetization below an apparent crit- 
ical temperature for unrealistically large systems. The 
relative magnetization persists below a finite tempera- 
ture for systems of iV = S^'"^ spins. For a practical im- 
plementation, the infinite system is a poor approxima- 
tion, and it remains poor even for finite systems that are 
astronomical in size. This conflicts with our notion of 
the thermodynamic limit, where the infinite system well 
describes crystalline solids of few billion unit cells JJJ] . 

First order finite-size scaling analysis is incomplete and 
fails to describe the slow decrease of T^^^ with the sys- 
tem size for all structures. For systems with no phase 
transition at finite temperatures, the shift away from Tc 
cannot be written as a simple power expansion in A^^/" . 
The partition function of the systems, and its first and 
second derivatives with respect to an external magnetic 
field, do not present critical points. They are continu- 
ous, well behaved, and show spontaneous magnetization 
for a broad range of temperatures. In the broad region 
near Tc=0, scaling properties of the magnetization and 
magnetic susceptibility satisfy power-laws as a function 
of N. 

The memory stabilizers presented in this work do not 
show a phase transition in the thermodynamic sense. 
However, for a wide range of temperatures and finite 
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size, there are many long paths of correlated spins that 
go through the structure resulting in a net macroscopic 
magnetization. These structures have free energy func- 
tions that spontaneously break global symmetry in the 
presence of environment-induced fluctuations, thus they 
stabilize the memory. 

We arrive to similar conclusions as in our previous 
work is]. Fig. [S] suggests that one way to increase the ap- 
parent critical temperature for a system of a given finite- 
size is by adding generators to each spin site. Struc- 
ture 3 is the best self-correcting memory as it has the 
broadest range of temperatures in which the system re- 
mains magnetized. The four simple two-body-interaction 
structures investigated have different levels of connectiv- 
ity. Wc find that the relationship between coordination 
number and the apparent finite size critical temperature 
Tmax is obvious. The number of generators is less 
important than the structure. The canonical stabilizer 
Hamiltonian seems to be thermodynamically a less sta- 
ble memory than the simpler pairwise based construction 
with the minimum number of generators (Structure 1), 
but it is more stable than the Ising chain. Kinetically, 
the canonical stabilizer could be the the most impervi- 



ous to fluctuations. For systems with all of the spins 
aligned, the lowest excited state energy for Structure 1 
is 2 J from the ground state, but this gap grows as 2fcJ 
for the canonical stabilizer. The multi-body interactions 
of the canonical stabilizer result in many large kinetic 
barriers that may be advantageous for preserving certain 
spin configurations. 

Finally, the exploration of stabilizer Hamiltonians de- 
fined by geometries of non-integer dimensions could yield 
self-correcting quantum memories with few multi-qubit 
interactions. Small finite size systems show an unusual 
order preservation for a broad range of temperatures 
making them suitable for a practical implementation of 
passive error correction. 
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Appendix A: Trees and trees with cycles 

For an Ising tree, there is only one path between any 
two nodes n and I. We define P/,„ as the product of 
Z operators of all of the nodes on the path between n 
and I on the spanning tree (including n and I). The 
transformation of Eq. [TT] can then be written succinctly 
as 

C/Z[„^^„];7t ^P^^^ (Ai) 
and M transforms to 

n n 

The expected value must be zero by symmetry and this is 
easy to confirm since each term contains the root factor 
Z[o 0] . Applying the transformation to the M operator 
of Eq. [51 results in 



n 



y ] ^[nd-i,n]^l^d-2,nd-i] 



n 



(A3) 
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Note that none of the terms contain ^[0,0] ^-nd the number 
of Z factors in each term is the distance between n and 
the root. 

The expectation vahie of the magnetization is the prod- 
uct of the polarization of ah the spins on the path from 
n to n\ in the free-spin basis. For the structures studied 
in Section IIII Al the spins that are in sibhng pairs have 
polarization a and otherwise have polarization e (except 
the root). This leads the relative magnetization to be 



(A4) 



where c„ is the number of spins that are in a pair between 
the root and n. One can express this graphically by la- 
beling every edge in the graph with an e if it is not part 
of a triangle, and with an a if it is part of a triangle (see 
Fig. [2]). One then starts from a node and multiplies the 
label of the edges between the node and the root on the 
spanning tree. Summing over all nodes yields Ea. ll2[ and 
a comparison to Eq. lA3l shows that Path(0, n) — {Pn,ni)- 
The magnetization operator squared, Af^ = M^, in 
the transformed basis is 

n.l 

= J2PojPo.n. (A5) 

The product of Pqj and Po,n results in the Z's that are 
in the intersection of the paths from the root to I, and 
from the root to n, to cancel. The Z operator on the 
last node in common is ^Last(i,n)j smd this node is in- 
cluded in the path from Z to n on the spanning tree. As 
an example, consider node 4 and node 9 of Fig. [51 Then 
PqA = ■2'[i,4]-^[o,i]-^[o,o], -fo,9 = Z[5,9]Z[i,5]Z[o,i]Z[ofi], and 
P4^g = 4]Z[]^ Q]Z[]^ 5]Z[5 9]. Node 1 is the last node 
in common and, as a result, /o,4-Po,9 = -P4,9-2^[o,i] f-i^d 
^[0,1] = ■2^Last(4,9)- These operators are the same for 
Fig. ^a) and (b), but the expectation values differ as 
explained below. 

Eg. lASl is simplified to: 



Last(/,n) ■ 



(A6) 



To calculate the expectation value of Pi^nZn,i, we must 
consider three cases. In the first case, the node n is 
contained in the path between the root and I or vice- 
versa, and ^Last(;,n) equals ^[„j,,„] or Z[i^j], respectively. 
The expectation value of Pi^nZn,i is simply the polar- 
ization of spins on the path from / to n excluding the 
node Last(Z,n). The polarization of each node is the la- 
bel of the edge connecting it to its parent and as a result 
{Pn,iZn,i) = Path(n, /). In the second case, the two nodes 
after the last node are not part of the same triangle. The 
polarization of spins at these nodes are independent and 
again {Pn,iZn,i} — Path(n, In the third case, the two 
nodes after the last node are part of the same triangle. 



The polarizations are not independent and two Z^s yield 
a single a. This is equivalent to taking a shortcut, and 
for all cases, (Pn.iZn^i) = Path(n, Z). Combining these 
observations with Eq. IA6I we obtain 



(A7) 



(M^) = ^Path(n,0. 



Appendix B: Canonical Stabilizer 

To define the full stabilizer, it is useful to start at the 
top level and work down. At level k there is one qubit 
(labeled 0), composed of three level k — 1 qubits labeled 
00, 01, and 02. We can define the logical operator as 



y{k) _ y(k-l) y(k-l) y(k-l) 

^0 ^ ^00 ^01 ^02 



and the highest order stabilizer elements as 

Aik) _ yi 
^01 — '^O 



_ (fe-l) (fc-l) 

.{k) _ (fe-l) (fe-l) 

^02 ^ ^00 ^02 



(Bl) 



(B2) 



Continuing this procedure, we define 



zif) 

V 




1 

^■ql 


-1)70-1) 


aU) 


_ yU-K 




-1) 


aU) 

^1)2 


— ^r,0 


^■q2 


-1) 



and one composite operator, 



(B3) 



(B4) 



where 7y is a fc — j + 1 string of trits. It is convenient 
to consider as a number in base 3. Also, we introduce 
four useful identities: 

(B5) 
(B6) 
(B7) 
(B8) 

Note how the products of A>^^ with Z^ interchange the 
1 and 2 labels. 

The Hamiltonian is 



40)70) 


— '^7)2 


-1) 


40)70) 




-1) 


40)70) 


- ^-qQ 


-1) 


4O) 4O) 


= /. 





fe 2 a'-'^^-i 

j—l x—l rj—O 



(B9) 



and there exists a unitary that transforms A's to single 
qubit Z's. The chosen unitary performs the following 
transformation: 

U-^rixU^ — Zj-ixou-'^y 

= ^))10U-i) -^7,200-1) I (BIO) 
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where is a string of I zeros and x = 1 or 2. Every 
physical qubit is denoted by a fc + 1 trit string with the 
first trit set to zero. This will transform the 3'^ — 1 stabi- 
lizer elements into S'^ — 1 single Z operators. To further 
specify the unitary, we set 

UZ^U^^Zok+i. (Bll) 

An explicit construction of U is as follows. Arrang- 
ing the spins on the tree that defines Structure 1 and 
starting at the spins at the maximum distance from the 
root, dmax, apply CNOT{nd^^^-i,nd^^J between ah 
spins connected on the graph at this distance. Then 
apply CNOT{nd^^^,nd^^^-i) to the same spins. Then 
move up one level and repeat the procedure first apply- 
ing CNOT{nd,„^^-2,nd^^^-i) at this distance and then 
CNOT{nd^^^^i,nd^^^-2)- Continue to the root. If we 
only used the CNOTs that are controlled by parents, 
we have the same unitary as Structure 1. Switching to 
CNOTs that are controlled by children allows us to trans- 
form logical Z's into Z's on single spins. The unitary is 
shown as a circuit for fc = 2 in Fig. |31 

W Zfj^U must be determined to calculate the magneti- 
zation. This is possible by constructing from stabi- 
lizer elements and Zq using Eqs. IB5[ IB61 and IB7I If we 
write /i = fikfJ-k-i-'-fJ-ifJ-o^ one can show that 

k 

Z, = 4n4.M.-....P.^.' (B12) 

where 

p.j = {2fij) mod 3 (B13) 
which exchanges I's for 2's and follows from Egs. lBSI and 

in 

Zq is the product of Z on all spins, and the element 
ZqAq^ is the product of all spins on a single k — 1 level 
qubit. Geometrically, each splits off one block of 3-'"^ 
spins from a block of 3'' level. Each A reduces the total 
number of spin operators by 1/3 until we reach a single 
qubit operator. 

To calculate the magnetic properties, we examine the 
transformed product of two Z operators, and Z^. As- 
sume that /i and v agree in the first q+l trits, the product 
is then 

k — q 

— ^/ifc/jfc_i.../ij_i^l/fci^fc_l...Pj_l (B14) 

After the transformation U Zi^iZy,U\ each remaining 
^r({i 2} corresponds to a single Z operator and each A^^q 
corresponds to two Z operators. The Z's will be unique 
between and v except for the only case where ry agrees. 



and The number of in- 

dependent Z's in U Z^Zi.U'^ denotes an effective distance 
between /i and 

5{n, v) = 2{k- q-l) + zeros(^fc_q_2.../xo) 
-t-zeros(z^fe_g_2---t'o) 

in which q -1- 1 is the number of leading trits that agree in 
/i and z^, and zeros(/ij.../io) counts the number of zeros 
in the string fij...fj,Q. The first trits that do not agree, 
^k-q-i and Vk-q-i, result in 1 or 2 independent Z's de- 
pending on whether either trit equals zero. For a given 
q, the minimal effective distance between the spins is 
2fc-2g- 1. 

Calculating the relative magnetization requires com- 
puting the effective distances between the root Zpt+i and 
all other spins Z^. Clearly, the number of leading trits 
that agree are the number of leading zeros in fi, L{fi). 
The result is 

(5(0^-+i, /x) = 2(fc - L(^) - 1) + (fc - - 1) 

+zeros(^fc_i(^)_2---/"o) + 1 
= 3fc - 4:L{fi) + zcros(/i) - 2. (B16) 

This effective distance serves the same role as the dis- 
tance from the root to the nodes in Structure 1. 

Fig. [7] shows the effective distance for each spin of the 
canonical stabilizer up to level k — 2. Each qubit at level 
k consists of three qubits of level fc — 1. The effective 
distances for each level k include effective distances from 
the central fc — 1 level qubit, and two equivalent sets 
of distances corresponding to the appended fc — 1 level 
qubits. The minimum effective distance in each new set 
is 3fc — 2 at the corner qubit, which occurs when there is 
exactly one zero in the string /i. In a new set, the effective 
distance increases by one for each additional zero in fi. 
Suppose /X labels a qubit in an appended block, then 
/ifc = and /Ltfc_i ^ 0. If ^ij = for j < fc — 1, then 
the qubit is effectively one step farther from the root 
than if /i^ = 1 or 2. Analogous to the Ising trees, the 
expected value of the magnetization is calculated using 
the polarization e and the effective distance of each qubit. 

With the minimum effective distance of 3fc — 2, the 
expected value of the magnetization is calculated to be 

(^M{k)) = (£{{k - 1)^ + 2(2 + ef-^e"^^-^, (B17) 

and the solution to this recursion formula is 

The 2 + e arises from the three choices of /ij and there 
are fc — 1 values of j in the appended block. 

The thermodynamic averag e of {M'^{k)) is obtained by 
performing a similar calculation but summing over the ef- 
fective distances between all spins. At level fc, each level 
fc — 1 qubit has the same internal magnetization squared. 
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FIG. 7. The effective distances of the canonical stabilizer up 
to fc = 2. For each k the black dot represents the root qubit. 
The corresponding is written next to each qubit. The 
numbers shown on the right are the effective distances from 
the root to the qubit at the same position. The bold box 
contains the effective distances inherited from the level k — 1. 



qubits in fc — 1 blocks, wc determine how many trits are 
required to distinguish a pair of physical qubits. It takes 
a single trit to determine which two logical blocks arc 
paired. Each block contains qubits labeled by fc — 1 dif- 
ferent trits. Similar to the magnetization, each trit con- 
tributes a factor of 2 -|- e, but now there are a total of 
2fc — 1 choices. The minimum effective distance between 
two sets is also 2fc— 1, since only the leading trit can agree 
(q=0). Overall, the thermodynamic average of (M^) is 
expressed as 



(M2(fc)) = 3(M2(fc-l)) 
+2(2 + e) 



(B19) 



and the solution to the recursion is given by 



fe-i 



(M2(fc)) = 3'= + 2.3*=-icE V 

■i-(CV3)'= 



= 3*^ 2 • 3*'"^C 



i-(CV3) 



(B20) 



(M'^{k — 1))- To compute the effective distances between 



where = {2 + e)e. 



